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Three possible scenarios of the statistical hadronization model are reexamined with the use of the 
Pt spectra of the PHENIX and very low pr PHOBOS measurements at -^/sjvjv = 200 GeV. These 
\£) ■ scenarios are: (a) full chemical non-equilibrium, (b) strangeness chemical non-equilibrium and (c) 

chemical equilibrium. Fits to the spectra are done within the Cracow single-freeze-out model, which 
takes into account both the expansion and resonance decays. Predictions for spectra of <f>, K(892)*° 
and 7T° are also given. The global variables like the transverse energy at midrapidity, the charged 
' particle multiplicity at midrapidity and the total multiplicity of charged particles are evaluated and 

| their predicted values agree qualitatively well with the experimental data. The thorough analysis 

/""N , within this model suggests that the chemical full non-equilibrium case is the least likely and both 

other cases are of the similar likelihood. It is also shown that if the full chemical non-equilibrium 
t-H ' freeze-out took place it could manifest itself in the enhancement of the 7r°-production at very low 

t-H , transverse momenta. 

fSJ . PACS numbers: 25.75.-q, 25.75.Dw, 24.10.Pa, 24.10.Jv 
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Q\ • I. INTRODUCTION 

o : 

Since the first run of the Relativistic Heavy Ion Collider (RHIC) the great amount of data on the hadron production 
from the hot and dense fireball created in a collision have been available. In this paper the application of the statistical 
i-^ | hadronization model (SHM) to the description of the fireball bulk properties [l[ is reexamined with the use of the pt 
I , spectra measured by the PHENIX Collaboration at the top RHIC energy y / ijvw = 200 GeV [|J . 

q ' In the SHM the formation process of each particle is described on the basis of the assumption that the accessible 
phase space is fully saturated (maximized) . Then the particle yields are determined by their phase space weight which 
is given by a statistical distribution (for a comprehensive review of the model see [1|). The main feature of this model 
is that it allows to deviate from the usually presumed chemical equilibrium of the fireball at the freeze-out. This 
has been achieved via the introduction of some new parameters, so-called phase-space occupancy factors: j q for light 
' quarks and 7 S for strange quarks in hadrons. In Ref. [l[ three possible cases were considered: 



1. Full chemical non-equilibrium, j q ^ 1, <y s ^ 1. 

2. Strangeness chemical non-equilibrium (semi-equilibrium), j q = 1, <y s ^ 1. 

3. Chemical equilibrium, j q = 1, <y s = 1. 

The phase-space occupancy factors 7 9 and 7 S together with the temperature T and baryon number chemical 
potential \ib comprise the full set of independent statistical parameters of the model. For all three cases of chemical 
non-equilibrium/equilibrium, values of these parameters have been determined in Ref. This was done for each 
centrality bin of the PHENIX measurement at ^/snn — 200 GeV from fits to the PHENIX identified hadron yields 
complemented with K* (892) / K~ and (f>/K~ ratios measured by the STAR Collaboration [1, 0, @|- 

However the possible expansion of the fireball is invisible in particle yield ratios (a collective flow is able to modify 
only the momentum spectra of a measured particle, but not its multiplicity). Therefore, in addition to the studies of 
particle ratios the analysis of the pr spectra is necessary to gain some quantitative information about the flow. Such 
analysis was done for the chemical equilibrium case of the SHM (the third point of the list above) in Ref. pj ■ In the 
present paper the similar analysis of the pr spectra will be performed for both chemical non-equilibrium cases of the 
SHM (the first and the second point of the list above). 

To describe the flow at the stage of the freeze-out, the single-freeze-out model of Refs. is applied. The model 

succeeded in the accurate description of ratios and px spectra of particles measured at RHIC. The main postulate of 
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the model is the simultaneous occurrence of chemical and thermal freeze-outs, which means that the possible elastic 
interactions after the chemical freeze-out are neglected. The conditions for the freeze-out are expressed by values of 
two independent thermal parameters: T and ub- The second basic feature of the model is the complete treatment of 
resonance decays. This means that the final distribution of a given particle consists not only of the thermal part but 
also of contributions from all possible decays and cascades. Feeding from week decays is included as well. Since in the 
original formulation [l,^ [l(| this model corresponds to the chemical equilibrium case of the SHM, the generalization 
of the single- freeze-out model to chemical non-equilibrium cases of the SHM will be done in the present paper. 

The global variables like the transverse energy at midrapidity (dEx / dr)\ m id.) , the charged particle multiplicity at 
midrapidity {dN c h/dri\ m id) and the total multiplicity of charged particles (N c h) are also evaluated for both chemical 
non-equilibrium cases of the SHM for different centrality bins of the PHENIX measurements at ^Snn = 200 GeV [ll| . 
These three variables are independent observables, which means that they are measured independently of identified 
hadron spectroscopy. Since model fits were done to identified hadron data (particle yield ratios and pt spectra) 
and the global variables are calculable in the generalized single freeze-out model, it was natural to check whether 
their estimated values agree with the data. This has proven to be true within 10% accuracy. It should be stressed 
also here that the centrality independence of the total multiplicity of charged particles per participant pair has been 
reproduced. The evidence for such scaling of the total multiplicity was reported by the PHOBOS Collaboration 12]. 

In some sense this work could be understood as an additional test of the correctness of the determination of the 
statistical parameters of the SHM since these parameters enter primordial distributions of hadrons in the fireball at 
the freeze-out. Thus fits of geometric parameters of the generalized single-freeze-out model are done with the use of 
the values of the statistical parameters obtained earlier in Ref. 1] and treated as input here. The general conclusion 
is that the best quality fits to the px spectra of identified hadrons are obtained for the strangeness chemical non- 
equilibrium case of the SHM. And what is surprising, in the chemical equilibrium case the spectra seem to be fitted 
better then in the scenario with the full chemical non-equilibrium. Whenever in this paper the term "strangeness 
chemical non-equilibrium case" or "full chemical non-equilibrium case" is used it means the case with its values of 
statistical parameters taken from Ref. [l[ and listed in Table|H Sec. lIV Al Additionally, the spectra of <fi and i^(892)*° 
resonances are predicted. In this way spectra of each particle species whose yield was used in determination of the 
statistical parameters of the model [l| are calculated in here. Also the measurement of the low momentum 7r° is 
proposed as a test, which could help to ascertain whether the full chemical non-equilibrium could happen in the 
fireball at the freeze-out or could not. Namely, values of j q determined in Ref. [l[ cause that the predictions for 
low-pT 7T° are about 40% greater in this case then in semi-equilibrium or equilibrium cases. 



II. THE SINGLE-FREEZE-OUT MODEL AND ITS GENERALIZATION 



The main assumptions of the model are as follows: (a) the chemical and thermal freeze-outs take place simultane- 
ously, (6) all confirmed resonances up to a mass of 2 GeV from the Particle Data Tables [l3| are taken into account, 
(c) a freeze-out hypersurface is defined by the equation 



t = y t 2 — r 2 . — Ty — r 2 = const , (1) 
(d) the four- velocity of an element of the freeze-out hypersurface is proportional to its coordinate 



(e) the following parameterization of the hypersurface is chosen: 

t = t cosh a\\ cosh a±, r x = r sinhaj^ cos</>, r y = r sinh a± sin <fi, r z = r sinh an cosh a±, (3) 
where otu is the rapidity of the element, a\\ — tanh _1 (r 2 /t), and aj_ controls the transverse radius: 



P = y r l + r l = T Sinn a ±- < Prnax , (4) 

where the restriction on the transverse size has been introduced, so p m ax gives the maximal transverse extension 
of the gas in the central slice during the freeze-out. This means that two new parameters of the model have been 
introduced, i.e. t and p max , which are connected with the geometry of the freeze-out hypersurface. 



3 



From Eq. |l| one can see that the beginning of the freeze-out process starts at t^P = r and r = in the c.m.s., 
which is also the laboratory frame in the RHIC case. At this moment the volume of the gas can be estimated as 

Vg = 27rrp 2 max , (5) 

which is simply the volume of a tube with a length 2r and a radius p m ax (2r is the maximal possible extension of the 
gas in the longitudinal direction at tP Q ). In the central slice the freeze-out ceases at tP = ^/t 2 + p max anc l ^ takes 
place at p = p max . 

The transverse velocity in the central slice can be expressed as a function of the transverse radius 

X (p) = tanh a x = -f== ■ (6) 
V T + P 

The maximum value of called the maximum transverse-flow parameter (or the surface velocity) is given by 

pmax Pmax _ Pmax I T 

V t2 + Pmax V /l + (Pmax/r) 2 

The invariant distribution of the measured particles of species i has the form Q 

dN,, 



d 2 p T dy 



p»da^f,(p-u), (8) 



where da^ is the normal vector on a freeze-out hypersurface, p ■ u = p^u^ , is the four-velocity of a fluid clement 
and fi is the final momentum distribution of the particle in question. The final distribution means here that fj is the 
sum of primordial and simple and sequential decay contributions to the particle distribution (for details see fiol. fl4j). 

For the most general case of the chemical non-equilibrium the primordial momentum distribution of particle species 
i is given by 

^primordial (2-5^ ~f~ 1) 1 , * 

h - ~7nZi~v?. \ r Z .. i ' (?) 



{2-Khcf 7 -i exp {S^| 



where Ei — (m 2 + p 2 ) 1 / 2 and m,i, pi, Si and gi are the mass, chemical potential, spin and a statistical factor of 
species i respectively. The chemical potential p,i = BipbB + Sips + I\P-i 3 , where Bi 1 Si and /| are the baryon 
number, strangeness and the third component of the isospin of the particle species in question, whereas /x's are the 
corresponding overall chemical potentials. The strangeness chemical potential ps is determined from the requirement 
that the overall strangeness equals zero. The chemical potential related to the third component of the isospin, pi 3 , is 
derived from the constraint that the charge to the net baryon ratio in the final state is the same as in the colliding 
nuclei. It has turned out that pj 3 is neglig ible at RHIC (| p h \< 1 MeV G3), so it will be omitted in further 
considerations. The non-equilibrium factor 74 reads 

(N*+N>) (N' a +Nl) rifV . 

7i = 7g 7s , (10) 

where 7 9 ( s ) is the light (strange) quark phase space occupancy factor, N l q and N l s are the numbers of light and strange 
quarks in the «th hadron, and N l = and iV| are the numbers of the corresponding antiquarks in the same hadron. 
With the use of Eqs. @ and (j3]), the invariant distribution ([5]) takes the following form: 



dN, 



d 2 Pr dy 

where 



J da (p ■ u) fi(p ■ u) = t 3 J J sinha^c?(sinha:^) J d£ (p ■ u) fi(p ■ u) , (11) 



p ■ u — rriT cosh («|| — y) coshaj_ — px cos£ sinha^ . (12) 
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III. TRANSVERSE ENERGY AND CHARGED PARTICLE MULTIPLICITY 

The experimentally measured transverse energy is denned as 



E T = Ej ■ sin I 



(13) 



where 9i is the polar angle, Ei denotes Ei — mjy (mjv means the nucleon mass) for baryons, Ei + mjv for antibaryons 



and the total energy Ei for all other particles, and the sum is taken over all L emitted particles [11 
The pseudorapidity density of particle species i is given by 



dN 1= f d 2 dy dNj 
drj J drj (Ppx dy 



= d p T — 



dNi 



Ei d 2 p T dy ' 

Analogously, the transverse energy pseudorapidity density for the same species can be written as 



(14) 



dE 7 



d z p T E t - — -f- 



p T dy dNi 



d PT PT ~fT 



Ei dNi 



drj J r± 1 p dr\ d?pT dy J Ei d 2 pT dy 

For the quantities at midrapidity one has (in the c.m.s., which is the RHIC case) 



(15) 



dNj 
drj 



mid 



PT 



d p T 



dNi 



vtiT d 2 pT dy ' 



(16) 



dE T ., 



mid 



{ J d 2 p T PT J^r y , * = baryon 

I d 2 p T PT J^r y , i = antibaryon 



J d 2 p T Pt 



d.N, 



i = others 



d 2 p T dy ' ' 

The overall charged particle and transverse energy densities can be expressed as 



(17) 



dNch 
drj 



dNj 
dtj 



(18) 



dE T 
drj 



sr^ dET.i 



(19) 



where A and B (B C A) denote sets of species of finally detected particles, namely the set of charged particles 
B = {ir + , 7r~, K + , K~ , p, p}, whereas A also includes photons, K^, n and n j45[. 

The total multiplicity of particle species i can be also derived (for the more formal proof see [lol ]) 



N 



J d 2 p T dy d 2^ dy = J d 2 PT dy J p^da^ fc(p ■ u) = J da J d 2 p T dy (p ■ u) fc(p ■ 
= J da J ^ (p ■ u) fi(p ■ u) = J da n t (T, ji B ,ls,l q ) = m(T, (i B , r fs,'rq) J da, 



(20) 



for any expansion satisfying the condition da^ ~ u,, on a freeze-out hypersurface and if the local statistical parameters 
are constant on this hypersurface (in the present model both conditions are fulfilled) . Note that the density of particle 
species i, rii, includes thermal and decay contributions. In practise the rapidity of the fluid element an should not 
be unlimited but should have its maximal value a™ aa: . Otherwise, the hypersurface volume and the total charged 
particle multiplicity would be infinite. Then one can express the hypersurface volume as 
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+ Q || Pmax/r 2tv 

Jda = r 3 J dct\\ J smha± d(8wha±) J d£ = 2w aJ ax Tp 2 max . (21) 

_ Q n,a J; 

Finally, the total multiplicity of charged particles can be obtained: 

N ch = 2tt af ax rp 2 max J2ni(T,^ B , 7snq ) = 2tt of ax rp 2 max n ch {T,^ B , ls , lq ) . (22) 

For a|p aa: the following reasonable assumption has been made: it is equal to the rapidity of leading baryons after 
the collision. This means that the fluid which has been created in the central rapidity region (CRR) could not move 
faster in the longitudinal direction then fragments of a target or a projectile after the collision. Therefore a™ ax should 
depend on the centrality of the collision, since the more central the collision is, the higher degree of the stopping of 
the initial baryons ought to happen in principle. There are two limiting cases, the maximum stopping happens for 
the most central collision whereas if the centrality approaches 100% the stopping disappears. Assuming additionally 
that u™ ax is a linear function of the centrality, the following parametrization can be derived (for details see Ref. [7[): 



af ax (c)=y p --^-(l-c), (23) 

where y p is the projectile rapidity, (Sy) the average rapidity loss and c is a fractional number representing the middle 
of a given centrality bin, i.e. c = 0.025 for the — 5% centrality bin, c = 0.075 for the 5 — 10% centrality bin, etc.. 
The BRAHMS Collaboration reports (Sy) = 2.05 for the 5% most central collisions at ^/s^n — 200 GeV (y p — 5.36) 



IV. RESULTS 



A. Determination of geometric parameters 



The determination of parameters of the model proceeds in two steps. First, statistical parameters T, /is, lq and 
7 S are fitted with the use of the experimental ratios of hadron multiplicities at midrapidity. This has been already 
done in Ref. [l| for all available centrality bins of the PHENIX measurements at ^/snn — 200 GeV 0. Having put 
values of these parameters into the theoretical expression for the invariant distribution, Eqs. (1111) and (|12p . two left 
parameters p m ax and r can be determined from the simultaneous fit to the transverse-momentum spectra of 7T , K , 
p and p. The fits are performed with the help of the \ 2 method. 

The final results for the geometric parameters p ma x and r are gathered in Table [J together with the corresponding 
values of x 2 /NDF for each centrality class additionally characterized by the number of participants N part . The results 
are given for all three cases of the SHM listed in Sec. Q] (for comparison the results for the chemical equilibrium case 
are repeated from Ref. Q). Other physical quantities like the surface velocity /3™ aa: , the volume at the beginning 

of the freeze-out vj 1 ^ and the maximal freeze-out time at the central slice tf a are also given there. Values of p m ax 

and r (therefore also vj^ and tj ) obtained in the case of full chemical non-equilibrium are substantially lower than 
corresponding values in the both other cases. This is because j s and j q are significantly greater than 1 in this case, 
so primordial densities given by Eqs. ([9]) and (fT0|) are also greater than in both other cases. And since fits are done 
to the same spectra, to keep the normalization unchanged, values of the geometric parameters have to decrease. 

Except the last three rows of Table U all fits have been done with the use of the pr spectra of identified charged 
hadrons measured by the PHENIX Collaboration in ^/sjvjv = 200 GeV Au-Au collisions @, [l7|. Centrality classes 
with footnote marks denote two bins for which fitted spectra are taken from Ref. [13]. These are 0—10% and 10 — 40% 
centrality bins and they arc not included in Ref. 0] , so values of the statistical parameters have not been fitted for 
them in Ref. [l|. But for these bins <j> meson spectra have been reported in Ref. [17[. Thus to make predictions for <p 
spectra, values of the statistical parameters have been taken as the averages of the values fitted for bins which added 
percent coverage equals — 10% or 10 — 40%. In the last three rows of TableQ]there are results of fits to the PHENIX 
data complemented with the low-pr data for tt^ 1 extracted from the PHOBOS measurements of (n + [Hj]. Since 

the particle ratio of 7t _ /tt + w 1 independently of pr and centrality (see e.g. Ref. ||), the low-pr values of n + and n~ 
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TABLE I: Values of the geometric parameters of the model for various centrality bins fitted with the use of the PHENIX final 
data for the pr spectra of identified charged hadrons 01, NDF=124. For bins with footnotes a and b the data are from Ref. [17| ]. 
Values of the statistical parameters are taken from pj. The last three rows show the results of fits to the set of data which 
include the PHENIX data and low-p T tt+ and vr" data taken as a half of the PHOBOS data for (ir + + tt") 0, here NDF=132. 



Centrality 


Npart 


T 


1MB 


7s 


7<J 


Pmax 


r 


jomax 
P± 




+ (2) 
l f.o. 


X 2 / NDF 


[%] 




[MeV] 


[MeV] 






[fm] 


[fm] 




[fm 3 ] 


[fm] 




0-5 


351.4 


141.1 


25.67 


2.430 


1.613 


7.24±0.09 


6.61±0.06 


0.74 


2177.5 


9.8 


0.74 


5-10 


299.0 


141.4 


24.52 


2.367 


1.61169 


6.82±0.08 


6.17±0.06 


0.74 


1804.4 


9.2 


0.73 


0-10" 


325.2 


141.25 


25.095 


2.3985 


1.6125 


7.03±0.08 


6.39±0.06 


0.74 


1985.8 


9.5 


0.80 


10-15 


253.9 


141.6 


25.27 


2.270 


1.603 


6.46±0.08 


5.81±0.06 


0.74 


1525.6 


8.7 


0.72 


15-20 


215.3 


140.8 


25.05 


2.266 


1.61497 


6.16±0.08 


5.48±0.05 


0.75 


1304.8 


8.2 


0.85 


20-30 


166.6 


141.0 


26.01 


2.212 


1.61387 


5.58±0.07 


4.96±0.05 


0.75 


969.7 


7.5 


1.24 


30-40 


114.2 


142.0 


25.75 


2.096 


1.608 


4.76±0.07 


4.33±0.05 


0.74 


617.0 


6.4 


1.64 


10-40 6 


171.8 


141.35 


25.52 


2.211 


1.61 


5.80±0.07 


5.17±0.05 


0.75 


1093.7 


7.8 


1.28 


40-50 


74.4 


141.7 


26.14 


2.003 


1.605 


4.06±0.06 


3.80±0.04 


0.73 


392.6 


5.6 


2.02 


50-60 


45.5 


141.0 


24.05 


1.876 


1.613 


3.39±0.06 


3.32±0.04 


0.71 


239.1 


4.7 


1.96 


60-70 


25.7 


140.2 


25.32 


1.636 


1.618 


2.72±0.05 


2.86±0.04 


0.69 


133.0 


3.9 


2.16 


70-80 


13.4 


141.7 


24.24 


1.026 


1.299 


2.40±0.06 


2.78±0.05 


0.65 


100.6 


3.7 


1.43 


0-5 


351.4 


154.6 


25.04 


1.231 


1.0 


8.35±0.10 


8.57±0.08 


0.70 


3752.1 


12.0 


0.57 


5-10 


299.0 


155.2 


24.73 


1.186 


1.0 


7.84±0.10 


7.97±0.08 


0.70 


3077.1 


11.2 


0.43 


0-10" 


325.2 


154.9 


24.885 


1.2085 


1.0 


8.10±0.09 


8.27±0.07 


0.70 


3407.3 


11.6 


0.56 


10-15 


253.9 


155.5 


26.29 


1.169 


1.0 


7.37±0.10 


7.41±0.07 


0.70 


2527.0 


10.4 


0.36 


15-20 


215.3 


154.6 


25.68 


1.147 


1.0 


7.07±0.10 


7.04±0.07 


0.71 


2212.2 


10.0 


0.39 


20-30 


166.6 


155.2 


27.18 


1.121 


1.0 


6.37±0.09 


6.32±0.06 


0.71 


1609.9 


9.0 


0.53 


30-40 


114.2 


155.7 


27.21 


1.080 


1.0 


5.47±0.08 


5.51±0.06 


0.70 


1036.6 


7.8 


0.78 


10-40 6 


171.8 


155.25 


26.59 


1.1293 


1.0 


6.65±0.08 


6.60±0.06 


0.71 


1833.5 


9.4 


0.53 


40-50 


74.4 


155.5 


26.74 


1.018 


1.0 


4.65±0.07 


4.82±0.06 


0.69 


654.6 


6.7 


1.07 


50-60 


45.5 


152.6 


21.62 


0.8906 


1.0 


4.06±0.07 


4.39±0.06 


0.68 


455.2 


6.0 


1.04 


60-70 


25.7 


152.2 


26.12 


0.8076 


1.0 


3.22±0.07 


3.73±0.05 


0.65 


243.0 


4.9 


1.32 


70-80 


13.4 


148.6 


23.82 


0.7163 


1.0 


2.59±0.06 


3.17±0.06 


0.63 


133.6 


4.1 


1.20 


80-92 


6.3 


150.8 


28.00 


0.6788 


1.0 


1.92±0.06 


2.69±0.06 


0.58 


62.3 


3.3 


1.21 


0-5 


351.4 


155.2 


26.4 


1.0 


1.0 


8.46±0.10 


8.84±0.08 


0.69 


3973.4 


12.2 


0.80 


5-10 


299.0 


155.2 


26.4 


1.0 


1.0 


7.99±0.10 


8.23±0.08 


0.70 


3302.6 


11.5 


0.61 


0-10" 


325.2 


155.2 


26.4 


1.0 


1.0 


8.23±0.09 


8.54±0.07 


0.69 


3629.8 


11.9 


0.80 


10-15 


253.9 


155.2 


26.4 


1.0 


1.0 


7.54±0.10 


7.67±0.08 


0.70 


2736.2 


10.8 


0.48 


15-20 


215.3 


155.2 


26.4 


1.0 


1.0 


7.11±0.10 


7.17±0.07 


0.70 


2275.5 


10.1 


0.48 


20-30 


166.6 


155.2 


26.4 


1.0 


1.0 


6.45±0.09 


6.47±0.07 


0.71 


1689.5 


9.1 


0.58 


30-40 


114.2 


155.2 


26.4 


1.0 


1.0 


5.57±0.08 


5.63±0.06 


0.70 


1097.2 


7.9 


0.77 


10-40 6 


171.8 


155.2 


26.4 


1.0 


1.0 


6.74±0.08 


6.76±0.06 


0.71 


1932.3 


9.6 


0.64 


40-50 


74.4 


155.2 


26.4 


1.0 


1.0 


4.68±0.07 


4.85±0.06 


0.69 


669.0 


6.7 


1.05 


50-60 


45.5 


155.2 


26.4 


1.0 


1.0 


3.83±0.07 


4.16±0.05 


0.68 


383.9 


5.7 


1.13 


60-70 


25.7 


155.2 


26.4 


1.0 


1.0 


2.99±0.06 


3.47±0.05 


0.65 


194.3 


4.6 


1.41 


70-80 


13.4 


155.2 


26.4 


1.0 


1.0 


2.22±0.06 


2.78±0.05 


0.62 


86.3 


3.6 


1.55 


80-92 


6.3 


155.2 


26.4 


1.0 


1.0 


1.71±0.06 


2.40±0.05 


0.58 


44.2 


2.9 


1.40 


0-15 c 


303.0 


141.4 


25.15 


2.356 


1.609 


6.82±0.08 


6.09±0.05 


0.75 


1778.0 


9.1 


0.87 


0-15 c 


303.0 


155.1 


25.35 


1.195 


1.0 


7.87±0.10 


8.00±0.07 


0.70 


3111.7 


11.2 


0.41 


0-15 c 


303.0 


155.2 


26.4 


1.0 


1.0 


8.02±0.10 


8.26±0.08 


0.70 


3337.0 


11.5 


0.59 



a Here statistical parameters are the averages of the parameters listed in two sequential rows above this row. 
'Here statistical parameters are the averages of the parameters listed in four sequential rows above this row. 

c Here statistical parameters are the averages of the parameters given for the — 5%, 5 — 10% and 10 — 15% centrality classes in the same 
case of the SHM. 
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FIG. 1: Values of the geometric parameters of the model from the seventh and eighth column of Table|T]for for the full chemical 
non-equilibrium case (7 a 7^ 1, j q 7^ 1). The lines are the best power approximations. 

spectra have been taken as one half of (ir + + ir~) reported by PHOBOS. However, some modification of the original 
PHENIX data 0] has been done to match the PHOBOS data conditions. Namely, the PHOBOS measurements were 
done for the 15% most central collisions (N part = 303), whereas the PHENIX ones for the — 5%, 5 — 10% and 
10 — 15% centrality bins. Since the treatment of counts includes the averaging over the number of events in a given 
centrality bin and for the same run the number of events in the 15% most central bin should be equal to the sum of 
numbers of events in the — 5%, 5 — 10% and 10 — 15% centrality bins, the rough approximation of the hypothetic 
measurement done in the — 15% centrality bin would be the average of the measurements done in the — 5%, 5 — 10% 
and 10 — 15% centrality bins. Such averages have been taken as the PHENIX data for the — 15% centrality bin. 
Also values of the statistical parameters taken for this case are the appropriate averages of the values given for the 
- 5%, 5 - 10% and 10 - 15% centrality bins. 

As it can be seen from the last column of Table U the best quality fits have been obtained for the strangeness 
chemical non-equilibrium case of the SHM. Also fits done in the chemical equilibrium case are slightly better than 
those presented for the full chemical non-equilibrium. This conclusion can be expressed in an informal quantifiable 
way by calculating the average of x 2 /NDF for each case of the SHM. So, for the chemical full non-equilibrium 
(x 2 /NDF) = 1.30, for the strangeness chemical non-equilibrium (x 2 /NDF) = 0.77 and for the chemical equilibrium 
(x 2 /NDF) = 0.9. Also the wider range of centrality fulfils the condition of the statistical significance, i.e. x 2 /NDF< 1, 
in the both cases of y q = 1 (up to 40% of centrality) than in the case of 7 9 7^ 1 (up to 20% of centrality) . Fits done 
with the inclusion of the low-pr i" measured by PHOBOS have confirmed the above conclusion, as it can be seen in 
three last rows of Table [J 

Values of the geometric parameters p max and r from Table U are presented in Figs.[T][2] as functions of N part . Also 
there the lines of the best power approximations are depicted, 

x ~ N parU x = Pmax, T, (24) 

with a scaling exponent n w 0.36 for p max and k ~ 0.28 for r. 

B. Identified hadron spectra 

Having obtained parameters of the model the spectra can be given with the use of Eqs. ([TT]) and (fT5|) . In Figs. [3] and 
|4] (top plots) the spectra of sums of negative and positive identified hadrons are depicted. This way of presentation 
is chosen to confront the model predictions for low-p-r values of spectra with the PHOBOS experimental data [l8[. 
Since the PHOBOS data are for the — 15% centrality bin, the PHENIX data for this bin have been simulated in the 
same way as explained in Sec. IIV Al 



8 



PHENIX 200 GeV 




p max =0.99774N part u 



*---- T=1.49392N parl 



100 



200 



300 



400 



N 



part 



FIG. 2: Values of the geometric parameters of the model from the seventh and eighth column of Table U for the strangeness 
chemical non-equilibrium case (7 S 7^ 1, y q = 1). The lines are the best power approximations. 




FIG. 3: The top plot presents invariant yields as a function of pr for RHIC at -^/sjvjv = 200 GeV. The PHOBOS data are for 
the 15% most central collisions with the error bars expressed as the sum of the systematic and statistical uncertainties 18]. The 
corresponding PHENIX data 0] are presented as the averages of the invariant yields for — 5%, 5— 10% and 10— 15% centrality 
bins. For the PHENIX data errors are about 10% and are of the size of symbols. Lines are the appropriate predictions of 
the single-freeze-out model for the full chemical non-equilibrium case (fit to the PHENIX data only). The bottom plot shows 
a deviation of data to the model, (fexp — ftheo)/crexp, where f e x P (theo) is the experimental (theoretical) value of the invariant 
yield at given pr and o cxp is the error of fexp- Both PHENIX and PHOBOS data are denoted by the same symbol for the 
same species, i.e. triangles are for (tt + + tt~), squares for (K + + K~) and open crosses for (p + p). 
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FIG. 5: Same as the bottom plot of Fig. [3] but for the chemical equilibrium case. 



In the case of chemical full non-equilibrium, Fig. [31 the low-p^ pions are mostly overestimated (w 33%), after then 
kaons (s» 21%) and the best predictions have been made for protons and antiprotons 5% above the data). In the 
case of chemical strangeness non-equilibrium, Fig.[H it is opposite, pions are predicted exactly, but kaons and protons 
and antiprotons are overestimated roughly equally (as 25%). For the chemical equilibrium case (see Fig. 4 in Ref. 0], 
it looks almost the same as the top plot of Fig. [?] here) the situation in the low-px range is similar to this in the 
chemical strangeness non-equilibrium case, namely pions are in complete agreement with the data, kaons are « 13% 
above and protons and antiprotons are the most overestimated, ~ 34%. The above discussion confirms the conclusion 
drawn from the comparison of the values of % 2 /NDF - the chemical strangeness non-equilibrium case seems to work 
in the best way as far as fits to the spectra are considered. 

To visualize the quality of fits and overall predictions, in the bottom plots of Figs. [3] and [4] and in Fig. [5] deviations 
of data to the model are presented. The deviation is defined as 



Jexp Jtheo /nr\ 

— . (25) 

where f eX p(theo) is the experimental (theoretical) value of the invariant yield at given px and a exp is the error of 
Jexp- The number of points which are entirely outside of the ±1 band is the greatest in the case of chemical full 
non-equilibrium (21, the most of them corresponds to pions with all the 1ow-j»t sample counted). In the both other 
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cases this number is the same and equals 8. Also the widest deviation is in the first case, it reaches -2.9, whereas in 
the case of chemical strangeness non-equilibrium the farthest point is -2.1 and in the chemical equilibrium case -2.8. 

To investigate into this problem from the other side, fits for the 15% most central bin have been done with the 
inclusion of the low-p^ 7r ± taken from the PHOBOS data fl8| . The data give values of (h + + h~) spectra of identified 
hadrons (h — ir, K, p) at very low px- But only for pions the particle ratio of h~ /h + 1 independently of px and 
centrality at RHIC (see e.g. Ref. 0). Thus values of tt^ 1 spectra for very low px are taken as one half of (n + + n~) 
reported by PHOBOS [Hj]. The results of fits have been gathered in the last three rows of Table HI The corresponding 
spectra are presented in the top plots of Figs. [B][H] In the bottom plots of Figs. [5][5] deviations of data to the model 
are depicted. Deviation figures shows explicitly what is expressed by the values of % 2 /NDF given in Table U- in the 
chemical strangeness non-equilibrium case spectra are fitted much better then in both other cases. The number of 
points outside the ±1 band equals 9 in this case and the farthest one is at -2.1 (see the bottom plot of Fig. [7]). In the 
chemical full non-equilibrium case this number is 19 and the farthest point is at -2.5 (see the bottom plot of Fig. 
In the chemical equilibrium case the number of points outside the ±1 band is 8, but one of them reaches the value 
-2.8 (see the bottom plot of Fig. [5]). 

From what has been explained so far one can see that the chemical full non-equilibrium freeze-out seems to be 
less likely in comparison with semi-equilibrium and equilibrium cases. And if j q = 1 indeed, both last cases will be 
practically undistinguishable, however semi-equilibrium will be in favor. 



C. 4> and if (892)*° spectra 

In this subsection the predictions for the spectra of <f> and if (892)*° resonances will be discussed. This is an 
interesting point since the yields of these resonances measured by the STAR Collaboration [J, 0, Q were used (with 
the basic yields of the identified hadrons measured by the PHENIX Collaboration 2]) to fit the statistical parameters 
of the model [l[ . In the fitting procedure presented here (to obtain the geometric parameters of the model p max and 
r, see Sec. IIV Ap . identified hadron spectra measured by PHENIX 0] have been explored. So the main source of the 
data used to test the SHM here and in Ref. [1] is the PHENIX measurement at ^/snn = 200 GeV. But predictions 
of the model should be compared with both PHENIX and STAR data, since the STAR data on K*(892)/K~ and 
4>/K~ ratios were also used in fits of the statistical parameters in Ref. [l|. However, one should keep in mind, when 





FIG. 8: Same as Fig. [7] but for the chemical equilibrium case. 
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FIG. 9: The top plot presents invariant yields of <f> meson measured via the K + K decay channel as a function of pr for the 
- 10% centrality bin at ^/sWn = 200 GeV. Data are from Refs. 0] (STAR) and [T| (PHENIX). The bottom plot shows a 
deviation of the data to the model predictions based on fits to the PHENIX spectra: chemical full non-equlibrium (circles), 
strangeness non-equilibrium (triangles) and chemical equilibrium (asterisks). 



the </> spectra are discussed, that the ~ 70% difference has been found between cf> yields at midrapidity measured by 
STAR pi] and PHENIX [17] for one common centrality been, — 10%. The reason for such behavior is still unknown 
and this is probably not a statistical fluctuation of the lowest mj PHENIX point, as suggested in Ref. [l[, because 
when the same mr range for both the STAR and the PHENIX <j> data is considered the difference still persists (see 
Ref. [59] in Ref. [Q). 

In Figs . 151 and ITTJ1 predictions for the <j> production are presented for — 10% and 10 — 40% centrality classes of the 
PHENIX measurement for all three cases of the SHM analyzed here. Additionally, the results for the equilibrium case 
but such that the statistical and geometric parameters of the model are fitted to the STAR data only are also depicted. 
This is the case considered in Ref. 0]: the statistical parameters (T — 160.0 MeV, /is = 24.0 MeV) are fitted to the 
STAR particle yield ratios (l9| an d the geometric parameters to the pt spectra of identified hadrons delivered by the 
STAR Collaboration in Ref. 2(j. Again, since the STAR identified hadron spectra [2Cj are for different centrality 
classes than the STAR 0-spectra |(|, the values of geometric parameters for — 10% and 10 — 30% centrality bins 
explored by STAR in meson measurements are the averages of the values fitted in Ref. Q for bins which added 
percent coverage equals — 10% and 10 — 30% respectively. This gives p ma x = 8.81 fm, r = 6.98 fm for the — 10% 
centrality bin and p ma x — 7.035 fm, r = 6.095 fm for the 10 — 30% centrality bin. Results corresponding to these two 
equilibrium (STAR) cases arc presented as long-dashed lines in Figs. 151 and [TUl Also both the PHENIX and the STAR 
data are depicted in these figures. Note that the STAR second bin is 10 — 30%, whereas the second PHENIX bin is 
10 — 40% (see Fig. [TU|) . Generally, as one can see in Figs. [51 and [TU1 all three cases of the SHM agree qualitatively 
with both the PHENIX and the STAR data when the predictions are based on the fits to the PHENIX spectra (solid, 
short-dashed and dashed lines). When the predictions are based on the fit to the STAR spectra (long-dashed lines) 
they agree with the STAR data only up to the intermediate transverse momentum range and overestimate the high 
Pt data. It is not clear why this happens, but within the PHENIX experiment the picture is consistent. 

Since Figs. [9l andfTOl are in the logarithmic scale, to see quantitative agreement or disagreement of each SHM case, 
one should go somehow to the linear scale. This has been done in the bottom plots of Figs . 151 and [TTJ1 where deviations 
of both the PHENIX and the STAR data to the model are depicted for all three cases of the SHM considered here. 
It is clearly seen that predictions in the equilibrium case are substantially better than in the both non-equilibrium 
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cases. In fact, the both non-equilibrium cases are practically ruled out, however the strangeness non-equilibrium case 
seems to behave slightly better than the full non-equilibrium case. 

In Fig. CCD results for (K*° + K*°)/2 spectra are presented together with the STAR data for - 10% and 10 - 30% 
centrality classes [|[ . This figure is very instructive since it explicitly shows that the data from different collaborations 
should not be mixed in any fitting procedure. As one can see predictions based on fitting to the PHENIX data for the 
- 10% centrality class (cf. Table [Q are very similar to each other and miss the STAR data on the [K*° + K*°)/2 
production mainly because of the different slope. In fact, these predictions are above the level of the STAR data for 
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Pt [GeV] 

FIG. 12: The top plot presents invariant yields of ty° in comparison with the half of (V + + n~) yields for RHIC at y / sjvjv = 200 
GeV. The PHOBOS data (diamonds) are for the 15% most central collisions [l||. The corresponding PHENIX data @] 
(triangles) are presented as the averages of the invariant yields for the — 5%, 5 — 10% and 10 — 15% centrality bins. The 
PHENIX data for tt° (stars) are for the 10% most central bin [24}. For both PHENIX data errors are about 10% and are of 
the size of symbols. Lines are the appropriate predictions of the single-freeze-out model. The bottom plot shows a ratio of 
predicted tv° to the half of measured (tt + +7r~) for the full chemical non-equilibrium case (triangles), the strangeness chemical 
non-equilibrium case (dots) and the chemical equilibrium case (open stars). Also values of the ratio calculated with the use of 
the experimental data are depicted (crosses). 

the — 10% centrality class in the low transverse momenta but they are below even the STAR data for the 10 — 30% 
centrality class in the high transverse momenta. Additionally the results for the equilibrium case but such that the 
statistical and geometricparameters of the model are fitted to the STAR data only are also depicted. This is the 
case considered in Ref. 0: the statistical parameters (T = 160.0 MeV, fi B = 24.0 MeV) are fitted to the STAR 
particle yield ratios [19j and the geometric parameters to the pr spectra of identified hadrons delivered by the STAR 
Collaboration in Ref. |20l|. Again, since the STAR identified hadron spectra [2(| are for different centrality classes 
than the STAR (K*° + K*°)/2 spectra @, the values of geometric parameters for - 10% and 10 - 30% centrality 
bins explored by STAR in (K* a + K*°)/2 measurements are the averages of the values fitted in Ref. Q for bins which 
added percent coverage equals — 10% and 10 — 30% respectively. This gives p max — 8.81 fm, r = 6.98 fm for the 
— 10% centrality bin and p max — 7.035 fm, r = 6.095 fm for the 10 — 30% centrality bin. Results corresponding 
to these two equilibrium (STAR) cases are presented as solid and shortest-dashed lines in Figs. [TTJ In fact some 
overestimation in normalization can be seen, mostly in the case of the — 10% centrality bin, but slopes are correct. 

D. 7T° spectra 

The occupancy factor "f q , when differs from one, could influence the 7r° spectra strongly. This is because for 7r° 
N q = Ng = 1 and N s = N s = 0. Then in the primordial distribution of ir° one has (see Eq. ©) 



(26) 
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FIG. 13: dNch/drj per pair of participants versus N par t for RHIC at y / ijvjv = 200 GeV. The original PHENIX data are from 
Ref. |IH. whereas the recalculated PHENIX data are from summing up the integrated charged hadron yields delivered in 
Ref. [21 • The lines connect the results and are a guide. 

where the chemical potential of pions is defined as 

^ = 2T\n lq . (27) 
It has turned out that in the case of chemical full non-equilibrium 



[i n ss ra T o (but [i n < m^a always) (28) 

for the fitted values of j q and T taken from Ref. [l| and listed in Table |U This means that in this case the values of 
the statistical parameters happen to hit the critical values for the Bose-Einstein condensation of neutral pions. This 
has been already stated by the authors of Ref. [l[ in Refs. [12, [23[ , namely that if 7 9 is freed from 1 but is kept in the 
range [l,7q r = e m T°/ 2T ], it goes to its critical value ^ during fitting procedure. If this happened really, this could 
enhance the production of ir°'s with very low pt- 

Predictions for ir° spectra are presented in Fig. [12] for two non-equilibrium cases of the SHM. The ir° spectrum in 
the chemical equilibrium case is roughly the same as that for the chemical strangeness non-equilibrium case, so it is 
not depicted. However it is impossible to compare the predictions for the low px with the data since the appropriate 
data have not been available yet. Thus in Fig.rT2]the comparison is done with the half of (ir + +7r~) spectrum delivered 
by PHOBOS for the 15% most central bin [18| and with the corresponding spectrum compiled from the PHENIX 
data [2|]. As it is shown in the bottom plot of Fig. the experimental ratio of 2ir° /(it + + ir~) w 1 in the range of 
Pt common for 7r ± and tt° [2l[ PHENIX measurements at ^/snn = 200 GeV, i.e. for 1 < pt < 3 GeV. 

One can see from the top plot of Fig. [T^] that down to the pr ~ 0.2 GeV all three cases of the SHM predict roughly 
the same spectrum of 7r° (the curve for the chemical equilibrium case is not depicted because, in the logarithmic 
scale, it would exactly cover the curve for the chemical strangeness non-equilibrium case). The difference between 
predictions in the chemical full non-equilibrium case and predictions in both other cases arises at very low transverse 
momenta and is about 40%. This can be seen very clearly in the bottom plot of Fig. [TSJ where the ratio of predicted 
7T° to the half of measured (ir + +tt~) is depicted as a function of pr for all three cases of the SHM. The enhancement 
of neutral pions over one half of charged pions is w 80% in the case of chemical full non-equilibrium, whereas for both 
other cases is w 30%. This suggests that the measurement of very low pt 7r°'s could be helpful to judge whether 
7<7 ~ lq r ( as i s claimed in Refs. [H, [23| on the basis of fits to particle yields/ratios) or 7 9 = 1. 



E. Transverse energy and charged particle multiplicity estimations 



The results of numerical estimations of dN c h/drj\ m id divided by the number of participant pairs for various centrality 
classes are presented in Fig. [T3] for RHIC at ^fsjf^ — 200 GeV. The results are given for two non-equilibrium cases of 
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FIG. 14: dEr/d-q per pair of participants versus N part for RHIC at ^snn = 200 GeV. The PHENIX data are from Ref. [ill ]. 
The lines connect the results and are a guide. 



the SHM (the estimates in the case of chemical equilibrium are almost the same as in the chemical strangeness non- 
equilibrium case and have been already presented in Ref. 0, see Fig. 6 therein). Additionally to the straightforward 
PHENIX measurements of the charged particle multiplicity density, the data from the summing up of the integrated 
charged hadron yields [1J are depicted in these figures, too (these data are called "recalculated", for more explanations 
see Ref. 0). Note that the recalculated data differ from the direct ones, especially for more central bins. This has 
been already noticed by the PHENIX Collaboration (see backup slides of [24]]). InJS the feeding of p(p) from A(A) 
decays is excluded. To diminish this effect, integrated p and p yields delivered in [2J were corrected to include back 
the feeding. The correction was done by the division by a factor 0.65, which is the rough average of a p^-dependent 
multiplier used by PHENIX Collaboration (see Fig. 4 in and Eq.(5) therein). 

Generally, in both presented cases the model predictions underestimate the directly measured (more) as well as 
recalculated (less) dN c h/ 'dr}\ m id- However, the estimates in the chemical full non-equilibrium case are slightly closer 
to the data and for the four most central bins they agree entirely with the recalculated data points. In the full 
range of centrality the predictions agree with the recalculated data within errors in this case and almost agree within 
errors in the chemical strangeness non-equilibrium case. In principle, since the fits of the geometric parameters of 
the model have been done to the same pt spectra here, which were integrated to deliver charged hadron yields in 
Ref. 0, the predictions for dN ' c h / 'dr]\ m id should agree exactly with the recalculated data. However, the transverse 
momentum spectra are measured in limited ranges, so very important Iow-Pt regions are blank in Ref. • To obtain 
integrated yields some extrapolations below and above the measured ranges are used. In fact these extrapolations 
are only analytical fits, but contributions from regions covered by them account for about 25 — 40% of the integrated 
yields 25]. These extrapolations could differ from the distributions obtained in the framework of this model and this 
could be the main source of the discrepancy between the predictions and the recalculated data. So the question why 
the significant underestimation of the predicted dN c h/drj\ m id with respect to the directly measured charged particle 
multiplicity density occurs should be addressed to the experimentalists rather: why does the directly measured 
dN c h / dri\ m id differ substantially from the sum of the integrated hadron yields for central collisions? 

The values of dET/drj\ m id per pair of participants as a function of participant pairs are shown in Fig. [14] for 
y/sNN = 200 GeV. The quality of the model predictions is much better in this case then for dN c h/dr]\ m id, they agree 
with the data almost completely. Note that predictions in both presented cases arc practically the same and do not 
differ from the corresponding results in the chemical equilibrium case (see Fig.7 in Ref. 0). 

Values of the ratio (dET/dr]}/(dN c h/dr]} as a function of N part are presented in Fig. [15] Again, as for dN c h/dr]\ m id, 
the values predicted in the chemical full non-equilibrium case are slightly closer to the data, they agree with the 
recalculated data within errors. For the most central bins both sets of predictions agree with the recalculated data 
within errors. (dEr / drj) / (dN c h/ drj) estimates done within the chemical equilibrium case are practically the same as 
in the strangeness non-equilibrium case (see Fig. 10 in Ref. 0). As far as the comparison with the direct data [ll| 
is concerned, the position of model predictions is very regular and exactly resembles the configuration of the data in 
each case, the estimates are only shifted up about 10% as a whole. 

The last discussed global variable is the total multiplicity of charged particles N c h, which can be calculated with 
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FIG. 16: N c h per pair of participants versus N par t for RHIC at ^snn = 200 GeV. The PHOBOS data are from Ref. [12| and 
the pp/pp data point of the UA5 measurement is from Fig. 39. 5 in Ref. [l3( ] . The lines connect the results and are a guide. 



the use of Eqs. (|22[) and ([23)) . The results presented as the total charged-particle multiplicity per participating pair 
versus N part are gathered in Fig. 1161 The both sets of predictions exhibit almost ideal centrality independence within 
the range of the PHOBOS measurement, i.e. N par t ~ 60 — 360. Note that in the chemical full non-equilibrium case 
also normalization agrees almost exactly with the data. In the case of chemical strangeness non-equilibrium the 6% 
underestimation has resulted but the predictions still agree with the data within errors. For the chemical equilibrium 
the similar underestimation was obtained (6.4%, see Fig. 12 in Ref. 0)- 

The general conclusion which could be drawn from the above discussion is that predictions for global variables agree 
pretty well with the data in each case of the SHM. However, the chemical full non-equilibrium scenario works slightly 
better in this respect. 
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V. CONCLUSIONS 

The extensive analysis of the RHIC data on the particle production in Au-Au collisions at ^/snn = 200 GeV has 
been performed within three possible scenarios of the statistical hadronization model. The SHM explored here is the 
generalized version of the model of Ref. [l| . The generalization means the explicite inclusion of the fireball expansion 
in a way as proposed in the single-freeze-out model of Refs. pi [sL [Tojl . 

Generally no definite rejection of any of these scenarios could be done on the basis of this analysis, since different 
observables prefer different scenarios. However, the chemical full non-equilibrium case seems to be the least likely. 
This is evident from the studies of particle spectra of stable charged hadrons, both from the statistical significance 
point of view (values of x 2 /NDF in Table |T| and from the behavior of the deviation factor (bottom plots of Figs. [3][5]). 
The ^-spectrum test confirms the above conclusion (see Sec. IIV Cj) . On the opposite, the global variable test prefers 
the full chemical non-equilibrium scenario and does not distinguish between strangeness chemical non-equilibrium and 
chemical equilibrium cases (see Sec. IIV B"|) but the differences are not significant. The semi-equilibrium and equilibrium 
scenarios seem to be of the similar likelihood, even though the 0-spectrum analysis discredits the strangeness chemical 
non-equilibrium case. This is because the fits to spectra of identified charged hadrons seem to weigh most when the 
conclusion is to be drawn from the present studies. The px spectra of stable charged hadrons comprise the most 
numerous and of the highest quality samples of the experimental data. For each centrality class of the PHENIX 
measurement at ^snn = 200 GeV, the px spectra counts more than 120 points, whereas each sample of the discussed 
resonance spectra or of the global variable has about 10 points. Also the measurement of a stable charged hadron 
seems to be more accurate since such a hadron is measured directly while resonances can be measured only via their 
decay products. The last point which supports the relevance of the fits to identified hadron spectra is that all data used 
in the fitting procedure (here particle yields and px spectra) and the data which predictions are compared with, should 
originate from the same experiment, as it has been shown explicitly in Sec. IIV CI in the example of (K*° + K*°)/2 
spectra (Fig. [TTj). This is the case of the fits to identified hadron spectra, since the statistical parameters were fitted 
to the sample comprising 6 particle yields from PHENIX , K ± , p and p) and only 2 yield ratios from STAR 
(see Ref. [lj). And with these parameters entering the expression for the invariant distribution, Eq. (|11[) . px spectra 
of 7r ± , K^, p and p measured by PHENIX have been fitted to determine the geometric parameters of the model. 
Thus the main results presented in Table U have been obtained within practically one experiment, i.e. the PHENIX 
Collaboration. 

The above remark should be kept in mind when the SHM predictions for yields of other particles (other then used 
in the fitting procedure) are compared with the data. For instance, in Ref. predictions for (anti)hyperons were 
done on the basis of fits from Ref. d (that is fits done to the data set of which the main part is from PHENIX at 
y/SNN = 200 GeV). However the conclusion is drawn from the comparison with the STAR data at ^Jsnn = 130 
GeV. The conclusion is that the chemical full non-equilibrium case is in favor. But in the main figure of Ref. [26[ 
(Fig. 2 there), which has led to this conclusion, there arc no corresponding predictions in the strangeness chemical 
non-equilibrium case. It is only stated there that these predictions are in between the chemical full non-equilibrium 
and the chemical equilibrium cases. All above arguments suggest that the strangeness chemical non-equilibrium 
case as well as the chemical equilibrium case have not been discredited entirely in the context of the (anti)hyperon 
production. 

Also particle yield fluctuations has been proposed as a definite test of what scenario of the SHM is the most likely [27| 
but one has to await for the appropriate data to make a conclusion. This test distinguishes between (semi)cquilibrium 
and non-equilibrium scenarios. But what is interesting, values of the statistical parameters given there for the chemical 
full non-equilibrium case are again at the condition for the Bose-Einstein condensation of neutral pions, 7 g « r ) q r 1 
(these values are T = 140.0 MeV and -f q = 1.62 [27|, which gives fi n = 135.079 MeV, Eq. (jST)). so fi n > m^o but this is 
the matter of rounding off (28|, if one takes j q = 1.619 then /i w < m x o). In fact, as it is explained in Refs. [13, [H[, j q r 
is the upper limit of the allowed range of 7 g superimposed before the fitting procedure has started. So by definition 
7 9 <2fq T always (if the value of j q put in a table of Refs. [2,[22,[23| happens to exceed 7^ r this is the result of rounding 
up [281 ] . as in the above-mentioned example). But from the technical point of view, when j q slightly exceeds this limit 
the fitting procedure will still proceed, since ir° yield is not included in the set of yields and/or ratios to fit. So the 
true upper limit should be e m ^ ± ^ 2T rather, because exceeding this limit causes divergences in primordial densities of 
7T + and 7r~, yields of which are included usually in the set of data to fit. Thus the fitted values of 7 g [l|, [22j, |23j seem 
to be not trustworthy. All these facts put in question the idea of introducing the parameter j q into the model. But 
this supports the conclusion that the chemical full non-equilibrium freeze-out is the least likely. Anyway, if values of 
7 9 were at the critical point for the Bose-Einstein condensation of ir° (as it is claimed in Refs. 0, [22|, HH), then the 
significant ^-overproduction at 1ow-pt could happen with respect to one half of (ir + + ir~) (see Sec. IIV D|) . which 
seems to be checkable at least in principle. 

And the last remark is that the present analysis has been done within a particular hypersurface, as given by Eqs. ([I])- 
(|1]). Of course, the natural question is to what extend the results depend on the choice of a hypersurface. One of the 
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indirect arguments pro this hypersurface are the results of fits to the PHENIX spectra of tt , K ± , p and p done in 
Ref. within the very popular blast-wave model (29j. For those fits x 2 /NDF« 3 — 4, so from the statistical point 
of view such a hypothesis should be rejected. On the opposite, in the present work x 2 /NDF< 1 has been obtained 
for all central and mid-central bins (see Table HJ . Thus at least for these bins the hypothesis that the hypersurface 
has the form as given here can not be rejected. Also the discussed possible overproduction of low-pr 7r°'s seems not 
to depend very much on a form of the hypersurface chosen since this effect is the result of approaching the condition 
for the Bose-Einstein condensation of neutral pions. This causes the abrupt increase of the distribution function of 
7T°> ^primordial ^ _ ^ when px — » for all hypersurfaces which have a region with the negligible flow. 

To summarize, in the view of this analysis the chemical full non-equilibrium freeze-out seems to happen least likely 
during Au-Au collisions at ^/sWn — 200 GeV and both other cases are of the similar likelihood. To help to verify 
this conclusion, the low-pr ir° measurement is proposed since at low pt the ratio of 7r° over one half of (ir + + 7r~) 
distinguishes very clearly between j q w 7^ r and j q = 1. 
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